Pulse and hold strategy for switching current measurements 
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We investigate by theory and experiment, the Josephson junction switching current detector in an 
environment with frequency dependent damping. Analysis of the circuit's phase space show that a 
favorable topology for switching can be obtained with overdamped dynamics at high frequencies. A 
pulse-and-hold method is described, where a fast switch pulse brings the circuit close to an unstable 
point in the phase space when biased at the hold level. Experiments are performed on Cooper pair 
transistors and Quantronium circuits, which are overdamped at high frequencies with an on-chip 
RC shunt. For 20 us switch pulses the switching process is well described by thermal equilibrium 
escape, based on a generalization of Kramers formula to the case of frequency dependent damping. 
A capacitor bias method is used to create very rapid, 25 ns switch pulses, where it is observed that 
the switching process is not governed by thermal equilibrium noise. 



I. INTRODUCTION 

A classical non-linear dynamical system, when driven 
to a point of instability, will undergo a bifurcation, where 
the system evolves toward distinctly different final states. 
At bifurcation the system becomes very sensitive and 
the smallest fluctuation can determine the evolution of 
a massive system with huge potential energy. This prop- 
erty of infinite sensitivity at the point of instability can 
be used to amplify very weak signals, and has recently 
been the focus of investigation in the design of quantum 
detectors to readout the state of quantum bits (qubits) 
built from Josephson junction ( JJ) circuits. Here we ex- 
amine in experiment and theory a pulse and hold strat- 
egy for rapid switching of a JJ circuit which is quickly 
brought near a point of instability, pointing out several 
important properties for an ideal detector. We focus on 
switching in a circuit with overdamped phase dynamics 
at high frequencies, and underdamped at low frequencies. 
This HF-overdamped case is relevant to experiments on 
small capacitance JJs biased with typical measurement 
leads. 

Classical JJs have strongly non-linear electrodynamics 
and they have served as a model system in non-linear 
physics for the last 40 years. More recently it has been 
shown that JJ circuits with small capacitance can also 
exhibit quantum dynamics when properly measured at 
low enough temperatures. Experimental demonstration 
of the macroscopic quantum dynamics in these circuits 
has relied on efficient quantum measurement strategies, 
characterized by high single shot sensitivity with rapid 
reset time and low back action. Some of these measure- 
ment or detection methods are based on the switching of 
a JJ circuit from the zero voltage state to a finite voltage 
stateji^fi Other detection methods are based on a dis- 
persive technique, where a high frequency signal probes 
the phase dynamics of a qubit i^&^ i 8 ' 9 ' 10 ! 11 ] 12 These dis- 
persive methods have achieved the desired sensitivity 



at considerably higher speeds than the static switching 
methods, allowing individual quantum measurements to 
be made with much higher duty cycle. In particular, 
the dispersive methods have shown that it is possible 
to continuously monitor the qubitJ^ However, for both 
static switching and dispersive methods, the sensitivity 
of the technique is improved by exploiting the non-linear 
properties of the readout circuit in a pulse and hold mea- 
surement strategy. This improved sensitivity of the pulse 
and hold method is not surprising, because when prop- 
erly designed, the pulse and hold technique will exploit 
the infinite sensitivity of a non-linear system at the point 
of instability. 

The general idea of exploiting the infinite sensitivity 
at an instable point is a recurrent theme in applications 
of non-linear dynamics. The basic idea has been used 
since the early days of microwave engineering in the well- 
known parametric amplifier— which has infinite gain at 
the point of dynamical instability. The unstable point 
can be conveniently represented as a saddle point for the 
phase space trajectories of the non-linear dynamical sys- 
tem. In the pulse and hold measurement method an ini- 
tial fast pulse is used to quickly bring the system to the 
saddle point for a particular hold bias level. The hold 
level is chosen so that the phase space topology favors a 
rapid separation in to the two basins of attraction in the 
phase space. The initial pulse should be not so fast that 
it will cause excessive back action on the qubit, but not 
so slow that it's duration exceeds the relaxation time of 
the qubit. The length of the hold pulse is that which is 
required to achieve a signal to noise ratio necessary for 
unambiguous determination of the resulting basin of at- 
traction. In practice this length is set by the filters and 
amplifiers in the second stage of the quantum measure- 
ment system. 

In this paper we discuss pulse and hold detection in 
the context of switching from the zero voltage state to 
the finite voltage state of a J J. We give an overview of 
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such switching in JJs, focusing on the HF-overdamped 
case. Switching detectors with overdamped high fre- 
quency phase dynamics are different from all other qubit 
measurement strategies implemented thus far, where un- 
derdamped phase dynamics has been used. However, the 
HF-overdamped case is quite relevant to a large number 
of experiments which measure switching in low capaci- 
tance JJs with small critical currents ! 15 ' 16 ! 17 ' 18 We show 
that by making the damping at high frequencies large 
enough, a favorable phase space topology for switching 
can be achieved^ In this overdamped situation the ex- 
ternal phase can be treated classically and contributions 
of macroscopic quantum tunneling (MQT) to the switch- 
ing probability can be neglected, in contrast to the un- 
derdamped case*? Experimental results are shown where 
on-chip RC damping circuits are used to create an HF- 
overdamped environment. We observe that for longer 
pulses of duration 20 vis, the switching process is initi- 
ated by thermal fluctuations in the overdamped system 
and thermal equilibrium is achieved at the base temper- 
ature of the cryostat (25 mK). For short pulses of du- 
ration < 25 ns, the switching is unaffected by thermal 
fluctuations up to a temperature of 500 mK, and the 
width of the switching distribution at low temperatures is 
rather determined by random variations in the repeated 
switch pulse. Although the detector apparently had the 
speed and sensitivity required for making a quantum 
measurement, we were unfortunately unable to demon- 
strate quantum dynamics of the qubit due to problems 
with fluctuating background charges. 



II. PHASE SPACE PORTRAITS 

The non-linear dynamics of a DC-driven J J can be pic- 
torially represented in a phase space portrait. We begin 
by examining the phase space portraits of the resistive 
and capacitively shunted junction (RCSJ), which is the 
simplest model from which we can gain intuitive under- 
standing of the non-linear dynamics. The RCSJ model 
consists of an ideal Joscphson clement of critical current 
io biased at the current level /, which is shunted by the 
parallel plate capacitance of the tunnel junction, Cj and 
a resistor R, which models the damping at all frequen- 
cies (see fig. Ufa)). The circuit parameters define the 
two quantities to p = y/Ia/foCj called the plasma fre- 
quency, and the quality factor Q — lo p RC. The dynamics 
is classified as overdamped or underdamped for Q < 1 or 
Q > 1, respectively. Here ipo = h/2e is the reduced flux 
quantum. 

The circuit dynamics can be visualized by the motion 
of a particle of mass <PqC in a tilted washboard poten- 
tial U(S) — —Ej(i5 + cos S) subjected to the damping 
force <Pq/R, where the particle position corresponds to 
the phase difference S across the junction and the tilt 
i is the applied current normalized by the critical cur- 
rent, i = I/Iq (see fig^a)). Below the critical tilt i = 1 
the fictitious particle will stay in a local minimum of the 
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FIG. 1: (a) RCSJ model of a Josephson junction, (b) Sim- 
ple model of a junction embedded in an environment with 
frequency dependent impedance. 



washboard potential (marked A in fig[^a)) correspond- 
ing to the superconducting state where V = (6) = 0. 
Increasing the tilt of the potential to i > 1, where local 
minima no longer exist, the particle will start to accel- 
erate to a finite velocity V — (8) > determined by the 
damping. If the tilt is then decreased below i = 1, the 
particle for an underdamped J J {Q > 1) will keep on 
moving due to inertia. Further decreasing the tilt below 
the level i < i r , where loss per cycle from damping ex- 
ceeds gain due to inertia, the particle will be retrapped 
in a local minimum. In terms of the current-voltage char- 
acteristic, this corresponds to hysteresis, or a coexistence 
of two stable states, V = Q and V > for a bias fixed in 
the region i r < i < 1. For the overdamped RCSJ model 
the particle will always be trapped in a local minimum 
for i < 1 and freely evolving down the potential for i > 1 
and there is no coexistence of two stable states. 

A phase space portrai t 19 ' 21 ' 22 of the RCSJ model is 
shown in fig. HJb). This portrait shows trajectories that 
the particle would follow in the space of coordinate (5) 
versus velocity (S) for a few chosen initial conditions. 
The topology of the phase space portrait is characterized 
by several distinct features. Fix point attractors marked 
"A" in fig. [U correspond to the particle resting in a lo- 
cal minimum of the washboard potential, and the saddle 
points marked "S" corresponds to the particle resting in 
an unstable state at the top of the potential barrier (com- 
pare with fig.^a)). Two trajectories surrounding A and 
ending at S are the unstable trajectories which define 
the boundary of a basin of attraction: All initial condi- 
tions within this boundary will follow a trajectory leading 
to A. We call this the 0-basin of attraction. The thick 
line B is a stable limiting cycle, corresponding to a free- 
running state of the phase 5, where the circuit is under- 
going Josephson oscillations with frequency lu — V/tpo- 
All trajectories leading to the limiting cycle B start in 
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the 1-basin of attraction, which is the region outside the 
0-basins. The existence of two basins of attraction in the 
phase space topology, and in particular the clear separa- 
tion of all 0-basins by the 1-basin, make the underdamped 
RCSJ circuit {Q > 1) appropriate for a switching current 
detector, as we discuss below. For the overdamped RCSJ 
circuit (Q < 1) attractors A and B do not coexist for any 
fixed bias condition, and it therefore can not be used for 
a switching current detector. However, the RCSJ model 
is not always the most realistic model for the dynam- 
ics of J J circuits, as the damping in real experiments is 
usually frequency dependent, and in the case of small ca- 
pacitance JJs, this frequency dependence can very much 
change the character of the damping. 

At high frequencies of the order of the plasma fre- 
quency of the junction (20-100 GHz for Al/AlOx/Al tun- 
nel junctions) losses are typically due to radiation phe- 
nomena, where the leads to the junction act as a wave 
guide for the microwave radiation. If we model the leads 
as a transmission line, the high frequency impedance 
would correspond to a damping resistance of the order 
of free space impedance Z w Z /2n = 60 fL With the 
small capacitance of a typical JJ as used in present ex- 
periments, this damping inevitably leads to overdamped 
dynamics Q < 1. It should be noted that for small ca- 
pacitance JJs, underdamped phase dynamics is hard to 
achieve in practice as high impedance all the way up to 
the plasma frequency is desired, and this requires an en- 
gineering effort where the high impedance leads need to 
be constructed very close to the junction.— However, at 
lower frequencies (typically below « 10 MHz) the junc- 
tion will see an impedance corresponding to the bias re- 
sistor R at the top of the cryostat, which can be chosen 
large enough to give Q > 1. The simplest circuit which 
captures the frequency dependence described above, is 
a JJ shunted by a series combination of a resistor R2 
and a capacitor C2 in parallel with the resistor Ri as 
shown in fig. HJb). At high frequencies where C2 is es- 
sentially a short, the circuit is described by the high- 
frequency quality factor Qi = oj p R\\Cj, where R\\ is 
the parallel combination of R\ and i?2- At low frequen- 
cies where C2 effectively blocks, the quality factor reads 
Qo = ujpRiCj. This model has been studied previously 
by several authors i 15 i 19 ' 24 i 25 ' 26 Casting such a circuit in a 
more mathematical language, it can be described by the 
coupled differential equation o 19 ' 25 
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where v — Vc/R\Iq is the reduced voltage across C2 
and p = R\C j I R2C2 reflects the value of the transition 
frequency, being u> ~ I/R2C2, between high- and low- 
impedance regimes. 

Phase space portraits for such a circuit are shown in 
figs- W[c) and (d). Here the y-axis shows the voltage v 




FIG. 2: (Color online) (a) Tilted washboard model, (b) Un- 
derdamped circuit biased at i r < i < 1. (c) Circuit with fre- 
quency dependent damping and Qi > Qi c and with Qi < Qi c 
(d). 
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which is directly related to 6. The topology of this phase 
portrait is also characterized by the coexistence of fix- 
points A and the limiting cycle B (not shown). However, 
for the parameters of fig. HJc) [Qq — 1.85, Q\ = 0.036 
and i — 0.8. ), the 0-basins and the 1-basin are now 
separated by an unstable limiting cycle C which does 
not intersect a saddle point. An initial condition which 
is infinitesimally below or above C will eventually end 
up either in an attractor A, or on B respectively. In 
fig- IH C ) we also see that the boundaries of the 0-basins 
are directly touching one another as a consequence of the 
existence of C. Thus, it is possible to have a trajectory 
from one 0-basin to another 0-basin, without crossing the 
1- basin. 

This same HF-overdamped model can however produce 
a new topology by simply lowering the high-frequency 
quality factor Q\. As we increase the high frequency 
damping, the unstable limiting cycle C slowly approaches 
the saddle points S. For a critical value of Q\ = Qi c , C 
and S will touch and the phase-portrait suddenly changes 
its topology. Fig.[^d) shows the phase space portrait for 
Qo = 1-85, Qi = 0.0065 < Q lc and i = 0.8 where we 
can see that C disappears and adjacent 0-state basins 
are again separated by the 1-basin - a topology of the 
same form as the underdamped RCSJ model. 



III. THE SWITCHING CURRENT DETECTOR 

The transition from the 0-basin to the 1-basin, called 
switching, can be used as a very sensitive detector. The 
idea here is to choose a "hold" bias level and circuit pa- 
rameters where the phase space portrait has a favorable 
topology such as that shown in figs. [2b) and (d). A 
rapid "switch pulse" is applied to the circuit bringing 
the system from A to a point as close as possible to the 
unstable point S. Balanced at this unstable point, the 
circuit will be very sensitive to any external noise, or to 
the state of a qubit coupled to the circuit. The qubit 
state at the end of the switch pulse can be thought of as 
determining the initial condition, placing the fictitious 
phase particle on either side of the basin boundary, from 
which the particle will evolve to the respective attrac- 
tor. The speed and accuracy of the measurement will 
depend on how rapidly the particle evolves away from 
the unstable point S, far enough in to the 0-basin or 1- 
basin such that external noise can not drive the system 
to the other basin. From this discussion it is clear that a 
phase space portrait with the topology shown in fig.[2lc) 
is not favorable for a switching current detector. Here 
the switching corresponds to crossing the unstable cycle 
C. Initial conditions which are infinitesimally close to C 
will remain close to C over many cycles of the phase, and 
thus a small amount of noise can kick the system back 
and forth between the 0-basins and the 1-basin, leading 
to a longer measurement time and increased number of 
errors. 

The measurement time is that time which is required 



for the actual switching process to occur and must be 
shorter than the relaxation time of the qubit. In the 
ideal case the measurement time would be the same as 
the duration of the switch pulse. A much longer time may 
be required to actually determine which basin the system 
has chosen. This longer detection time is the duration of 
the hold level needed to reach a signal to noise ratio larger 
than 1, which is in practice determine by the bandwidth 
of the low-noise amplifier and filters in the second stage of 
the circuit. In our experiments described in the following 
sections, we used a low noise amplifier mounted at the 
top of the cryostat which has a very limited bandwidth 
and high input impedance. While this amplifier has very 
low back action on the qubit circuit (very low current 
noise), it's low bandwidth increases the detection time 
such that individual measurements can be acquired only 
at < 10 kHz repetition rate. Since many measurements 
(10 4 ) are required to get good statistics when measuring 
probabilities, the acquisition time window is some 0.5 
seconds and the low frequency noise (drift or 1// noise) 
in the biasing circuit will thus play a role in the detector 
accuracy. 



IV. FLUCTUATIONS 

The measurement time of a switching detector will de- 
pend on fluctuations or noise in the circuit. The phase 
space portraits display the dissipative trajectories of a 
dynamical system, but they do not contain any informa- 
tion about the fluctuations which necessarily accompany 
dissipation. For a switching current detector, we desire 
that these fluctuations be as small as possible, and there- 
fore the dissipative elements should be kept at as low a 
temperature as possible. Analyzing the switching current 
detector circuit with a thermal equilibrium model, we can 
calculate the rate of escape from the attractor A. This 
equilibrium escape rate however only sets an upper limit 
on the measurement time. When we apply the switch 
pulse, the goal is to bring the circuit out of equilibrium, 
and we desire that the sensitivity at the unstable point 
be large enough so that the measurement is made before 
equilibrium is achieved (i.e. before thermal fluctuations 
drive the switching process). 

Equilibrium fluctuations can cause a J J circuit to jump 
out of it's basin of attraction in a process know as ther- 
mal escape. The random force which gives rise to the es- 
cape trajectory will most likely take the system through 
the saddle point S, because such a trajectory would re- 
quire a minimum of energy from the noise source^ For 
the topology of phase space portraits shown in figs. [2jb) 
and (d), thermal escape will result in a switching from 
a 0-basin to the 1-basin, with negligible probability of a 
" retrapping" event bringing the system back from the 1- 
basin to a 0-basin. However, for the topology of fig. [UJc), 
thermal escape through the saddle point leads to another 
0-basin, and thus the particle is immediately retrapped 
in the next minimum of the washboard potential. This 
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process of successive escape and retrapping is know as 
phase diffusion, and it's signature is a non-zero DC volt- 
age across the JJ circuit when biased below the critical 
current, i < f . 

Phase diffusion can occur in the overdamped RCSJ 
model, or in the HF-overdamped model when parame- 
ters result in a phase space topology of fig. [2^c) . In the 
latter switching process can be identified which 

corresponds to the escape from a phase diffusive state to 
the free running state, or to crossing the unstable lim- 
iting cycle C in fig. |^c) which marks the boundary be- 
tween the phase diffusive region and the 1-basin. This 
basin boundary C is formed by the convergence of many 
trajectories leading to different S, and the escape pro- 
cess of crossing this boundary is fundamentally differ- 
ent than escape from a 0-basin to the f-basin. Numer- 
ical simulation o 19 ' 27 ' 28 of switching in JJs with such a 
phase space topology show that escape over the unsta- 
ble boundary C is characterized by late switching events, 
which arise because even a small amount of noise near 
this boundary can kick the system back and forth be- 
tween the f-basin and the many 0-basins for a long time 
before there is an actual escape leading to the limiting 
cycle B. 

The rate of thermal escape from a 0-basin can be cal- 
culated using Kramers' formula 22 ' 29 ! 30 

r = K ^exp(~A£;/fc B T), (2) 

with AE being the depth of the potential well from A to 
S, ks Boltzmann's constant and T the temperature. The 
prefactor kujo/2it is called the attempt frequency, where 
k < I is a factor which depends on the damping. An- 
alytical results for k were found by Kramers in the two 
limiting cases of underdamped (Q > I) and overdamped 
(Q < 1) dynamics. For the application of Kramers' es- 
cape theory we require that AE fc^T, i.e. thermal 
escape is rare, so that each escape event is from a ther- 
mal equilibrium situation. The fluctuations in thermal 
equilibrium are completely uncorrelated in time, which is 
to say that the strength of the fluctuations are frequency 
independent (white noise). Furthermore, the Kramers 
formula assumes absorbing boundary conditions, where 
the escape process which leads to a change of the basin 
of attraction has zero probability of return. These condi- 
tions restrict the direct application of Kramers formula in 
describing switching in JJ circuits 3 ^ to the case of the un- 
derdamped RCSJ model such as that depicted in fig.^b). 
In principle one could apply Kramers formula to the over- 
damped RCSJ model, where the escape is from one well 
to the next well (switching between adjacent attractors 
A), but experiments thus far are unable to measure a 
single 2tt jump of the phase, as this corresponds to an 
extremely small change in circuit energy. 

Thermal induced switching of small capacitance 
Josephson junctions which experience frequency depen- 
dent damping as modeled by the circuit of Fig. [3b), was 
analyzed in experiment and theory by the Quantronics 



group 15 ' 26 who generalized Kramers result. The theoreti- 
cal analysis was subject to the constraint that the dynam- 
ics of the voltage across the shunt capacitor v is under- 
damped (i.e. the quality factor a = R2R\\C'2lo/ i Po{Ri + 
R%) 3> 1 where R\\ is the parallel resistance of R\ and 
Ra) so that the dynamics of v is subject to the fast-time 
average effects of the fluctuating phase S. Separating 
timescales in this way, the switching of v could then be 
regarded as an escape out of a meta-potential, B, formed 
by the averaged fluctuating force in the tilted washboard 
potential F = i — (sin 5) — v. Assuming non-absorbing 
boundary conditions, this "escape over a dissipation bar- 
rier" can be written as a generalization of Kramers' for- 
mula 

Here D(v) is the position-dependent diffusion constant, 
and B = J^*(F/XD)dv, where Vf, and v t stand for the 
bottom and the top of the effective barrier, respectively. 
Detailed expressions can be found in refs i 15 ' 32 

In section I VII we use these escape rate formulas to 
analyze pulse and hold switching measurements. We 
demonstrate that long switching pulses lead to thermal 
equilibrium switching, whereas short pulses switch the 
circuit in a way that is independent of temperature at 
low temperatures, with the switching distribution deter- 
mined by noise in the switch pulse rather than noise from 
the cooled damping circuit. 

V. EXPERIMENTS 

Experiments investigating junction current-voltage 
characteristics (IVC) as well as pulsed switching behav- 
ior were carried out in a dilution refrigerator with 25 mK 
base temperature. A block diagram of the measurement 
setup is shown in fig. El^a). A low noise instrumentation 
amplifier (Burr-Brown INA110, noise temperature 1.3 K 
at 10 kHz) is measuring the voltage across the sample 
while the sample is biased by a room temperature volt- 
age source either via the bias capacitor Cb, or in series 
with a bias resistor Ri,. The capacitor bias method was 
used for experiments with fast current pulses of duration 
t p = 25 ns, while the conventional resistor bias method 
was used for long pulse experiments with r p = 20 us, as 
well as for IVC measurements. 

Three different samples are discussed in this paper 
which differ primarily in the range of the measured 
switching current (3 nA to 120 nA), and in the type 
of circuit used for the damping of the phase dynamics. 
These different damping circuits are labeled in the order 
in which they were implemented, and are represented in 
fig.^a) as the blocks Fi, F 2 , F 3 . These environments can 
be modeled as RC filters with different cut-off frequen- 
cies, as schematically be represented in fig. |3Jb). 

The key parameters for each sample are given in Table 
HI Sample I consisted of a Cooper pair transistor (CPT) 
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Sample 


Type 


CPT 


Shunt JJ 


Fi 


F 2 


F 3 


Ej/Ec 


la 


I SW 


Ej/Ec 


Jo 


1 ' sw 


Hi 


Ci 


R 2 


c 2 




c 3 


I 


CPT 


32.9 


58.5 


3 








60 


1 










II 


CPT 


29 


51.6 


4.2 








60 


1 


7.2 


0.24 






III 


Quantronium 


2.2 


21 


12 


30.3 


158 


120 


1000 


3 


7.2 


0.24 


600 


1.4 



Currents are in [nA], resistances in [fi] and capacitances in [nF] 

TABLE I: An overview over the parameters for the three different samples. Filter Fi resembles the cryostat leads or a cold 
SMD filter. Filter F2 is the on-chip damping circuit and F3 is an on-chip RC-filter. 
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FIG. 3: (a) Block diagram of the experimental setup, (b) 
Schematic diagram of the model used to describes the three 
different filters F\,F2,F^ which define the damping circuit, 
(c) Micrograph of the sample showing the two on-chip RC 
filters F\ and F2. (d) Magnified view of the center part of 
the chip, F2 and the bias capacitor, (e) Electron microscope 
picture of the Quantronium (sample III). 



embedded in an environment denned solely by the twisted 
pair leads of the cryostat which is modeled as F\ . Sam- 
ple II was a CPT fabricated in parallel with Sample I, 
having nearly identical parameters, differing only in that 
sample II was embedded in a micro-fabricated on-chip 
HF-damping circuit F2. Sample III is a Quantronium 1 
embedded in the same HF-damping F2 used with sam- 
ple II, but with an additional micro-fabricated on-chip 
low-pass filter F3 . The on-chip RC-environments F2 and 



F3 used for samples II and III and the bias capacitor 
were fabricated with a two-step optical lithography pro- 
cess. The capacitors were actually two capacitors in se- 
ries, formed by a plasma-oxidized Al ground plane cov- 
ered with a Au top plate. The top plates are connected 
to the rest of the circuit via resistors which are formed 
from the same Au film as the top plate, having a typ- 
ical sheet resistance of 1.2 VL/U. The capacitors of F3 
could be measured quite accurately, from which we ob- 
tain a specific capacitance of 13.6 fF /um 2 that is used to 
determine all on-chip capacitors. Figure [3jc) shows the 
essential parts of the chip and the components defining 
the high-frequency environment. The bright rectangular 
area on the left side is the top plate of the capacitor, 
and the thin leads leading to the right are the resistors of 
filter F3. Figure E^d) shows in detail the biasing capaci- 
tor Cb The bright trapezoidal area on the left is the top 
plate of the capacitor C2 and the areas surrounded by 
dashed lines are damping resistors R2/2. Figure [3je) is 
an electron microscope picture showing the A1/A1 2 03/A1 
tunnel junctions, which were fabricated in a third layer of 
electron beam lithography, with the standard two-angle 
evaporation through a shadow mask. Figure Ele) shows 
the quantronium circuit of sample III. 
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FIG. 4: IV curves of sample I without designed RC environ- 
ment (solid line) and of sample II with specially designed RC 
environment (dashed line). 
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The effect of the on-chip HF damping from F2 on the 
phase dynamics as can be seen in fig. 2] where the IVC 
of samples I (solid line) and sample II (dashed line) are 
shown. These two CPT samples differ essentially by the 
presence of F2 in sample II. We sec that the typical phase 
diffusion shape of the IVC— of sample I is absent in 
sample II which shows a sharp supercurrent and hys- 
teretic switching. The presence of the on-chip environ- 
ment in sample II effectively reduces phase diffusion as 
can be explained by a phase-space topology as shown in 
figure [Hd). However, the very low value of I sw is a di- 
rect indication of excessive noise in the circuit. Therefore 
the on-chip low-pass filter F3 was implemented in sample 
III, improving the switching current to a value of 75% of 
the critical current. In the remainder of this paper we 
concentrate on investigating the switching behavior of 
sample III. 

The ability to suppress phase diffusion opens up the 
possibility to study fast switching with HF-overdamped 
phase dynamics for the first time. We used the pulse and 
hold method to measure switching probabilities of sample 
III as a function of the amplitude of the switch pulse for 
two cases: A long pulse of 20 us where the switching was 
found to be controlled by equilibrium thermal escape, 
and a short pulse of 25ns, where the switching is clearly 
a non-equilibrium process. 

The long pulses were formed by applying a square volt- 
age pulse through the bias resistor. The response to a 
simple square pulse is shown in fig. EJa), where the ap- 
plied voltage pulse is shown, and several scope traces of 
the measured voltage over the CPT are overlayed. Here 
we see that the switching causes an increase in the volt- 
age over the sample which can occur at any time during 
the applied pulse. In order to do statistics we want to 
unambiguously count all switching events. Late switch- 
ing events are difficult to distinguish from non-switching 
events as the voltage does not have time to rise above the 
noise level. We can add a trailing hold level as shown in 
fig. H^b). This hold level and duration must be chosen 
so that there is zero probability of switching on the hold 
part of the pulse. The response to such a pulse shows 
that it is now easy to distinguish switch from non-switch 
events. In this case the hold level is used simply to quan- 
tize the output, and the switching which occurs during 
the initial switch pulse is found to be a thermal equilib- 
rium escape process as discussed below. 

The fast pulses were formed with a new technique 
where a voltage waveform consisting of a sharp step fol- 
lowed by linear voltage rise is programmed in to an arbi- 
trary waveform generator. The slope on the sharp step 
(dV / dt) pu i se is typically 6-7 times larger than the linear 
rise during the hold, (dV/dt)hoid- The voltage waveform 
is propagated to the chip through a coax cable having 
negligible dispersion for the sharp 25 ns voltage step used. 
An on-chip bias capacitor Cj will differentiate the voltage 
waveform to give a sharp current pulse followed by a hold 
level, I = CbdV/dt, which is shown in fig.[Sfc). From the 
measured step amplitude needed to switch the junction 



(a) (b) (c) 




20 40 60 80 20 40 60 80 25 50 75 100 
t [us] t Qis] t [|xs] 



FIG. 5: (a) Square pulse and response, (b) Switch pulse 
with hold level and response, (c) Switch pulse with hold level 
and response, generated by the capacitive bias method. The 
finite voltage in the non-switch case is due to the series filter 
resistance and the two-point measurement setup. 



and the value of C& = 1.4 pF, we calculate a pulse ampli- 
tude of 360 nA through C&. Due to the symmetry of the 
filter stages F\ to F3, only half of this 25 ns pulse cur- 
rent flows through the junction, with the other half flow- 
ing through the filter to ground. Thus the peak current 
through the junction during the 25 ns pulse I p — 180 nA, 
which is larger than Iq. Exceeding Iq for this very short 
time is not unreasonable, bearing in mind that the circuit 
is heavily overdamped at high frequencies, and a strong 
kick will be needed to overcome damping and bring the 
phase particle close to the saddle point. 

The hold level for these fast pulses is 40 ps,very much 
longer than the switch pulse, and its duration is set by 
the time needed for the response voltage to rise above 
the noise level. The rate of this voltage rise depends on 
the hold current level because after the switch we are es- 
sentially charging up the second stage filter and leads, F$ 
and Fi, with the hold current, Ihold = Cb{dV/dt)hold — 
56 nA. For the low level of hold current used in these ex- 
periments, we can follow the voltage rise at the junction 
with the 100 kHz bandwidth low noise amplifier at the 
top of the cryostat. Typically we turn off the hold cur- 
rent and reset the detector when the sample voltage is 30 
uV, so that the junction voltage is always well below the 
gap voltage V2A = 400 pV, and therefore quasi-particle 
dissipation during the hold can be neglected. 

Pulsed switching measurements were performed were 
a sequence of 10 3 to 10 4 identical pulse-hold-reset cycles 
was applied to the sample while recording the voltage re- 
sponse of the sample. A threshold level was used to dis- 
tinguish switching events (1) from non-switching events 
(0) as depicted in fig. [|3a) . The maximum response volt- 
age achieved during each cycle is found and a histogram 
of these values is plotted as seen in fig. \$[b). The hold 
level and duration are adjusted so as to achieve a bimodal 
distribution in the histogram, with zero events near the 
threshold level, meaning that there is zero ambiguity in 
determining a switch event from a non-switch event. We 
further check that the hold level itself, without the lead- 
ing switch pulse, gives no switches of the sample. The 
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sequence of switching events is stored as a binary se- 
quence Yi in temporal order. From this sequence we can 
calculate the switching probability, 



(a) 



P 



1 - 



and the auto-correlation coefficients, 



Y)(Y i+k -Y) 



Yf 



(5) 



where fc is the "lag" between pulses. The auto-correlation 
is a particularly important check for statistical inde- 
pendence of each switching event. A plot of for 
k = 1...1000 is shown in figure [He) and the randomness 
and low level of indicates that all switching events are 
not influenced by any external periodic signal. When the 
circuit is not working properly, pick up of spurious sig- 
nals up to the repetition frequency of the measurement, 
clearly shows up as a periodic modulation in the auto- 
correlation 7-fc. Of particular importance is the correla- 
tion coefficient for lag one r\ which tells how neighboring 
switching events influence one another. Fig. [DJd) shows 
r\ as a function of the wait time t w between the end 
of the hold level and the start of the next switch pulse. 
For large values of t w , n fluctuates around not exceed- 
ing 0.05, which shows that any influence of a switching 
or non-switching event on the following measurement, is 
statistically insignificant. As t w is decreased however, a 
positive correlation is observed, with r% increasing expo- 
nentially with shorter t w . Positive correlation indicates 
that a switching event (a "1") is more likely to be fol- 
lowed by another switching event. Fig. [B^d) shows a fit 
to correlation r\ to the function 



r\ = 3.345 • exp 



33.3ps 



(6) 



We can extrapolate the fit to the time r* = 40.25 ps 
where the auto-correlation becomes n = 1, meaning that 
once the circuit switches it will always stay in the 1-state. 
In our experience, increasing the capacitance of filter F\ 
causes r* to increase, from which we infer that the in- 
crease in the correlation r% for short t w is resulting from 
errors where the detector is not properly reset because it 
does not have time to discharge the environment capaci- 
tance before a new pulse is applied. For the experiment 
shown in figure [S] the time constant of the environment 
was estimated to be 3 ps. These observations indicate 
that it is necessary to bring the junction voltage very 
close to zero before the retrapping will occur, and the 
detector will reset. For good statistics many pulses are 
required and a short duty cycle is desirable in order to 
avoid effects from low frequency noise as discussed sec- 
tion III. By studying the correlation coefficient r\ in this 
way, we can choose an optimal duty cycle. 
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FIG. 6: (a) Response of the sample to a pulse sequence result- 
ing in switches (1) and no-switches (0) of the sample, (b) Peak 
voltage obtained during a current pulse, indicating good sep- 
aration between the switch and non-switch signal, (c) Auto- 
correlation function r*,. (d) Correlation coefficient n vs. the 
wait time, with fitted exponential decay function. 



VI. ANALYSIS 

The switching probabilities were thus measured and 
the dependence on the amplitude of the switch pulse, 
P(Ip) was studied as as a function of temperature. Each 
measurement of P(I P ) began with a pulse sequences hav- 
ing pulse amplitude resulting in a switching probabil- 
ity P — 0, and the pulse amplitude was successively 
increased until P = 1. The measurement produces an 
" S-curve" as shown in figure [JJ where the experimental 
data for the long pulse duration t p — 20 ps is shown with 
crosses. The S-curves were taken at temperatures 100, 
200, 300, 400 and 500 mK (right to left) respectively. 

We compare the measured data to theoretical predic- 
tions based on thermal escape as discussed in section IV. 
The filter Fi causes a rounding of the applied square 
voltage pulse, which is accounted for by calculating the 
escape probability for a time dependent current^, 



P = 1 - exp - 



1 



di/dt 



T(i')di' 



(7) 



where the escape rate T can be found using either eqns.[2] 
or [3] The simulated S-curves using eqn. [3] are plot- 
ted in figure [JJ as solid lines for the temperatures cor- 
responding to the measured data. Sample parameters 
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FIG. 7: Switching probability as a function of pulse height 
in the temperature range T = 100, 200, 300, 400, 500 mK 
(right to left) for a pulse duration t p = 20 us. Crosses show 
measured data and simulated data is shown as solid lines. 



used for this calculation are the measured bias (including 
filter) resistance R\ = 11600 f2, the measured high fre- 
quency damping resistor, R2 — 7.2 ft, the high frequency 
damping capacitance C2 = 0.207 nF, the junction capac- 
itance Cj — 30 fF, and the calculated critical current 
Iq = 148 nA. The critical current Iq = 148 nA is not 
the bare critical current Iq = 158 nA since the quantron- 
ium was biased at a magnetic field such that a persistent 
current of « 10 nA was flowing in the loop. These param- 
eters are all independently determined, and not adjusted 
to improve the fit. However, the capacitance of filter F\ 
was uncertain, having a nominal value of lOnF, and un- 
known temperature dependence below 4 K. Cooling the 
same capacitance to 4 K, we observed a decrease of C\ by 
around 10%. This capacitor C\ determines the rounding 
of the square voltage pulse, and thus the time depen- 
dence of the current applied to the junction. We found 
that it was necessary to assume C\ = 3 nF in order for 
the simulations to agree with experiment. This low value 
of C\ at low temperatures is not unreasonable, as circuit 
simulations with the nominal value of 10 nF showed that 
the initial pulse would not exceed the hold level, which 
clearly is not possible because excellent latching of the 
circuit was observed. 

From the experimental and simulated S-curves, we de- 
fine the switching current of the sample as the pulse am- 
plitude that gives 50 % switching P(I SW ) — 0.5 and the 
resolution is defined from the S-curve by AI = I p {P = 
0.9) — Ip{P = 0.1). A comparison of experimental and 
theoretical I sw vs. T and AI/I SW vs. T is shown in 
fig. [51 We see that the experimental data for the long 
pulses (points marked by an X) are in reasonably good 
agreement with the simulated values when the theory of 
switching in an environment with frequency dependent 
damping is used (escape from a meta-potential, equation 
[3]) which is plotted as a solid line in fig. [8ja). We note 
that for the 20 ps pulses, escape occurs at bias currents 
i ~ 0.7, where the phase space has a topology as shown 
in figure HJd). Hence we can neglect phase diffusion and 



escape is from a saddle point, so that the non- absorbing 
boundary condition assumed in the theory is valid. For 
comparison, we use the overdamped Kramers formula 
(equation [5]) to simulate the S-curve and calculate I sw 
and AI, which is shown by the dashed line in fig.[8j Here 
the prefactor k(Q) is given in reb2£ and we have used the 
high frequency quality factor Qi = 0.027 as determined 
by the resistor R2 only. We see that the Kramers for- 
mula overestimates I sw by some 25% (fig. [8ja)) and is 
worse than the simulation based on cqn. [31 in reproduc- 
ing the temperature dependence of AI (fig.[Sfb)). In fact, 
the experimental data for the 20 ps pulses only shows a 
weak increase in AI over the temperature range studied, 
whereas both theoretical curves predict a slight increase 
in AI. Thus an equilibrium thermal escape model ex- 
plains the data for long, 20 ps pulses reasonably well 
and the data is better explained by the theory of escape 
with frequency dependent damping, than by the simpler 
theory embodied in the overdamped Kramers formula. 
However the correspondence with the former theory is 
not perfect. We may explain these deviations as being 
due to the fact that the quality factor a = 4.49 (see sec- 
tion [TV]) does not really satisfy the condition for validity 
of the theory, a ~S> 1. 

Experimental data for the short pulses of duration 
T p = 25 ns generated by the capacitive bias method is 
plotted in fig. [5] as circles. Here we see that the value 
of I sw is constant in the temperature range studied, in- 
dicating that escape is not from a thermal equilibrium 
state. For the ideal phase space topology, as shown 
in figure (2^d), the initial pulse would bring the phase 
particle arbitrarily close to the saddle point S for the 
hold bias level. If the separation in to the basins of at- 
traction occurs before thermal equilibrium can be estab- 
lished, we would not expect temperature dependence of 
I sw . In this case, the width of the switching distribution 
will be determined not by thermal fluctuations, but by 
other sources of noise, such as random variation in the 
height of the switch pulse. These variations are signifi- 
cant because the 1/f noise from the waveform generator 
must be taken in to account when generating the train of 
pulses over the time window of the measurement which 
was about 0.5 sec. In our experiments however, we may 
not have achieved a constant hold level since the volt- 
age ramp from the waveform generator is not perfectly 
smooth. Knowing the bias capacitor we can calculate an 
average hold level of ihoid = 0.35, somewhat lower than 
the critical value of ihoid = 0.67 necessary to achieve the 
phase space topology of figure [2][d). Nevertheless, we 
observe excellent latching of the circuit for these 25 ns 
switch pulses. We conclude that the observed tempera- 
ture independence of I sw , and the fact that I sw exceeds 
In by 20% is consistent with a very rapid switching of the 
junction. 

We can rule out excessive thermal noise as a reason for 
the temperature independent value of I sw for the short 
pulses. By measuring the gate voltage dependence of I sw 
as a function of the temperature, a clear transition from 
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FIG. 8: Switching current normalized to critical current (a) 
and relative resolution (b) of sample III. Crosses indicate mea- 
sured values for r p = 20 us, solid and dashed lines are calcu- 
lated values using a generalization of Kramers' large friction 
result and Kramers' original result, respectively. Circles are 
measured values for t„ = 25 ns. 



2e to e periodicity was observed in the temperature range 
250 mK to 300 mK. For the size of the superconducting is- 
land used in this experiment, we can estimate a crossover 
temperature T* rj 300 mK, above which the free energy 
difference between even and odd parity goes to zero. 34 
Hence, we know that the sample is in equilibrium with 
the thermometer below T* , and therefore heating effects 
that might occur in the short pulse experiments, can not 
explain the fact that the observed I sw is independent of 
temperature. 

Thus we have achieved a very rapid, 25 ns measure- 
ment time of the switching current, which should be suffi- 
cient for measurement of the quantum state of a quantro- 
nium circuit. For qubit readout, not only the measure- 
ment time is important, but also the resolution of the 
detector. For the 25 ns pulse, we obtained the resolution 
of AI/I SW = 0.055, or AI = 9.9 nA. This implies that 
single shot readout is possible for a Quantronium with 
parameters Ec — 0.5 K and Ej/Ec — 2.5 where the 



switching current of the two qubit states at the optimal 
readout point differ by 9.6 nA. Numerous experiments 
were made with microwave pulses and continuous mi- 
crowave radiation to try and find the qubit resonance. 
However, due in part to uncertainty in the qubit circuit 
parameters (level separation) and in part to jumps in 
background charge, no qubit resonance was detected in 
these experiments. 

VII. CONCLUSION 

Fast and sensitive measurement of the switching cur- 
rent can be achieved with a pulse-and-hold measurement 
method, where an initial switch pulse brings the JJ cir- 
cuit close to an unstable point in the phase space of 
the circuit biased at the hold level. This technique ex- 
ploits the infinite sensitivity of a non-linear dynamical 
system at a point of bifurcation, a common theme in 
many successful JJ qubit detectors. We have shown that 
with properly designed frequency dependent damping, 
fast switching can be achieved even when the high fre- 
quency dynamics of the JJ circuit are overdamped. With 
an on-chip RC damping circuit, we have experimentally 
studied the thermal escape process in overdamped JJs. A 
capacitor bias method was used to create very rapid 25ns 
switch pulses. We demonstrated fast switching in such 
overdamped JJs for the first time, where the switching 
was not described by thermal equilibrium escape. The 
methods presented here are a simple and inexpensive way 
to perform sensitive switching current measurements in 
Josephson junction circuits. While we have shown that 
the sensitivity can be high, the effect of back-action of 
such a detector is still unclear and might be a reason 
why no quantum effects were observed. In contrast to 
the readout strategy presented here, all other working 
qubit-readout strategies, both static switching and dis- 
persive, are based on underdamped phase dynamics. 
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